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Abstract 

We compute the 1 = 2 vrvr scattering length at pion masses of ~ 294, 348 and 484 MeV in fuhy- 
dynamical lattice QCD using Liischer's finite-volume method. The calculation is performed with 
domain-wall valence-quark propagators on asqtad-improved MILC configurations with staggered 
sea quarks at a single lattice spacing, b ~ 0.125 fm. Chiral perturbation theory is used to perform 
the extrapolation of the scattering length from lattice quark masses down to the physical value, 
and we find 771,^02 = —0.0426 it 0.0006 it 0.0003 it 0.0018, in good agreement with experiment. 
The 1 = 2 TTTT scattering phase shift is calculated to be 5 = —43 it 10 it 5° at |p| ~ 544 MeV for 
~ 484 MeV. 
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I. INTRODUCTION 



The scattering of two pions at low energies is the simplest of all dynamical strong-interaction 
processes from the theoretical standpoint. As pions are the pseudo-Goldstone bosons asso- 
ciated with the spontaneous breaking of chiral symmetry, their low-momentum interactions 
are constrained by the approximate chiral symmetries of Quantum Chromodynamics (QCD), 
and the scattering lengths for vrvr — > vrvr in the s-wave are uniquely predicted at leading order 
in chiral perturbation theory (x-PT). Higher orders in the chiral expansion give contribu- 
tions to these scattering lengths that are suppressed by powers of (m^/A^)^, where m^r is 
the mass of the pion, and A^ ~ 1 GeV is the scale of chiral symmetry breaking. However, 
part of these higher-order contributions are from local counterterms with coefficients that 
are not constrained by chiral symmetry alone. Lattice QCD is the only known technique 
with which one can rigorously calculate these strong-interaction quantities. A lattice calcu- 
lation of the TTTT scattering length as a function of the light-quark masses, rrig, will allow the 
determination of the relevant counterterms that appear in the chiral lagrangian. Of course, 
this is only true for values of m^r within the chiral regime where -C A^. Even though 
it is likely that in the future lattice calculations will be performed at the physical values of 
the light-quark masses, and while a chiral extrapolation of observables will not be necessary, 
the values of the counterterms will be essential as the chiral lagrangian will still be used to 
compute more-complex strong-interaction processes that are too costly to compute directly 
on the lattice. A second important motivation for studying tttt scattering with lattice QCD 
is its impact upon weak-interaction processes (and physics beyond the Standard Model), 
such as K ^ 271, and the determination of fundamental constants of nature associated with 
electroweak physics, such as CP violation in the Standard Model. 

Lattice QCD calculations are performed on a Euclidean lattice, and the Maiani- Testa 
theorem demonstrates that S-matrix elements cannot be determined from n-point Green's 
functions computed on the lattice at infinite volume, except at kinematic thresholds 0|. 
However, Liischer has shown that by computing the energy levels of two-particle states in 
the finite- volume lattice, the 2 — 2 scattering amplitude can be recovered The energy 
levels of the two interacting particles are found to deviate from those of two non-interacting 
particles by an amount that depends on the scattering amplitude and varies inversely with 
the lattice spatial volume. A number of lattice QCD calculations of vrvr scattering in the 
1 = 2 channel have been performed previously using Liischer's method. For obvious reasons, 
calculations were initially performed in quenched QCD 0, 0, S IE H, S E 13, IE H 
H H H El H H El but in a recent tour-de-force study, the CP-PACS collaboration 



exploited the finite- volume strategy to study 1 = 2, s-wave vrvr scattering in fully-dynamical 
lattice QCD with two fiavors of improved Wilson fermions |2j], with pion masses in the 
range m^r ~ 0.5 — 1.1 GeV. 

A lattice QCD calculation of tttt scattering should fulfill several conditions: (i) quark 
masses small enough to lie within the chiral regime so that the extrapolation to the physical 
point is reliable, (ii) a lattice volume large enough to avoid finite-volume effects from pions 
"going around the world" but small enough for the energy shifts to be measurable and (iii) 
a lattice spacing small enough so that discretization effects are under control. Based on 
these requirements we have performed a fully-dynamical lattice QCD calculation with two 
degenerate light quarks and a strange quark (2+1) with pion masses in the chiral regime, 
~ 294 MeV (319 configurations), 348 MeV (649 configurations) and 484 MeV (453 
configurations). The configurations we have used are the publicly-available MILC lattices 
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with dynamical staggered fermions of spatial dimension L ~ 2.5 fm, and we have used 
domain-wall (isl . l24l l25l. propagators generated by LHPC ^ at the Thomas Jefferson 
National Laboratory (JLab). 

This paper is organized as follows. In SectionlHlwe discuss Liischer's finite- volume method 
for extracting hadron-hadron scattering parameters from energy levels calculated on the 
lattice. In Section lllll we describe the details of our mixed-action lattice QCD calculation. 
We also discuss the relevant correlation functions and outline our fitting procedures. In 
Section |TV| we present the results of our lattice calculation, and the analysis of the lattice 
data with x-PT. In Section El we conclude. 



II. FINITE- VOLUME CALCULATION OF SCATTERING AMPLITUDES 

The s-wave scattering amplitude for two particles below inelastic thresholds can be deter- 
mined using Liischer's method j|2j] , which entails a measurement of one or more energy levels 
of the two-particle system in a finite volume. For two particles of identical mass, m, in 
an s-wave, with zero total three momentum, and in a finite volume, the difference between 
the energy levels and those of two non-interacting particles can be related to the inverse 
scattering amplitude via the eigenvalue equation p 



(1) 



where 6{p) is the elastic-scattering phase shift, and the regulated three-dimensional sum is 



S(r^) = 



|j|<A 



AttA 



(2) 



The sum in eq. (j21) is over all triplets of integers j such that | j | < A and the limit A ^ oo 
is implicit jl^] . This definition is equivalent to the analytic continuation of zeta-functions 
presented by Liischer 0. In eq. ((T)), L is the length of the spatial dimension in a cubically- 
symmetric lattice. The energy eigenvalue En and its deviation from twice the rest mass of 
the particle, AEn, are related to the center-of-mass momentum p„, a solution of eq. (^, by 



AE„ 



E„ 



2m 



— 2m 



(3) 



In the absence of interactions between the particles, |pcot5| = cxo, and the energy levels 
occur at momenta p = 27rj/L, corresponding to single-particle modes in a cubic cavity. 
Expanding eq. (Q) about zero momenta, p ~ 0, one obtains the familiar relation ^ 



Ana 



a fa 
1 + C- + c.(- 



(4) 



^ We thank the MILC and LHP collaborations for very kindly allowing us to use their gauge configurations 

and light-quark propagators for this project. 
^ We have chosen to use the "particle physics" definition of the scattering length, as opposed to the "nuclear 

physics" definition, which is opposite in sign. 
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with 



1 'j'^"^ 1 1 1 
ci = - E - 4A = -2.8372 , ca = c? ^ E = 6.37518 , (5) 

and a is the scattering length, defined by 

a = hm^^ . (6) 

For the J = 2 tttt scattering length, 02, that we consider in this work, the difference between 
the exact solution to eq. ([Q) and the approximate solution in eq. (@]) is much less than 1%. 
However, in determining the phase-shift associated with the first excited state on the lattice, 
the full eigenvalue equation in eq. (H)) is solved. 



III. DETAILS OF THE CALCULATION 

A. Lattices, Actions and Propagators 

Our computation of the / = 2 vrvr scattering amplitude consists of a hybrid lattice QCD 
calculation using staggered sea quarks and domain- wall valence quarks [28,, 2^ lUl ^■ 
The parameters of the three sets of Nf = 2 + 1 asqtad-improved ji^, MILC configura- 
tions generated with staggered sea quarks that we used in our calculations are shown 
in Table HJ In the generation of the MILC configurations, the strange-quark mass was fixed 
near its physical value, brng = 0.050, (where b is the lattice spacing) determined by the 
mass of hadrons containing strange quarks. The two light quarks in the three sets of con- 
figurations are degenerate, with masses bmi = 0.007, 0.010 and 0.020. These lattices were 
HYP-blocked [371 138^ in order to avoid large residual chiral symmetry breaking. We used 
the domain-wall valence propagators that had been previously generated by the LHP col- 
laboration on each of these sets of lattices. The do main- wall height is m = 1.7 and extent 
of the extra dimension is L5 = 16. As this is a mixed-action calculation, the parameters 
used to generate the light-quark propagators have been "matched" to those used to gen- 
erate the MILC configurations. This was achieved by requiring that the mass of the pion 
computed with the domain-wall propagators be equal (to few-percent precision) to that of 
the hghtest staggered pion computed from staggered propagators generated with the same 
parameters as the given gauge configuration [36]. The parameters used in the generation 
of the domain-wall propagators are shown in Table |l] A Dirichlet boundary condition at 



TABLE I: The parameters of the MILC gauge configurations and LHPC domain-wall propagators 
used in om calculations. For each propagator the extent of the fifth dimension is L5 = 16. 



Config Set 


Dimensions 


brrii 


hrus 


bmdwf 


lO'^ X brrires 


# configs 


MILC_2064f21b679m007m050 


20^ X 64 


0.007 


0.05 


0.0081 


1.604 ±0.038 


319 


MILC_2064f21b679m010m050 


20^ X 64 


0.010 


0.05 


0.0138 


1.552 ±0.027 


649 


MILC_2064f21b679m020m050 


20^ X 64 


0.020 


0.05 


0.0313 


1.239 ± 0.028 


453 



^ Technical aspects of hybridization are discussed in Ref. ,32j and ,3^ . 
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the midpoint of the time direction {t = 32) has been used in generating the LHPC propa- 
gators in order to reduce the cost of nucleon matrix-element calculations without affecting 
their accuracy^. Unfortunately, this is not the case for light-pseudoscalar mesons where the 
signal is sustained for longer time intervals, and hence both the systematic errors and the 
statistical errors can be improved using correlators of longer time extent. 

B. Correlators, Projections and Fitting Methods 

In order to perform our calculations of the vrvr correlation functions we used the programs 
Chroma and QDP++ written at JLab under the auspices of SciDAC [s^. In this program- 
ming environment, a few lines of C++ code were required to construct the two distinct 
propagator contractions that contribute to vrvr interactions in the 1 = 2 channel. As it is 
the difference in the energy between two interacting pions and two non-interacting pions 
that provides the scattering amplitude, we computed both the one-pion correlation func- 
tion, CT,(t), and the two-pion correlation function C^^(p, t), where t denotes the number of 
time-slices between the hadronic-sink and the hadronic-source, and p denotes the magnitude 
of the (equal and opposite) momentum of each pion. 
The single-pion correlation function is 

= J2(^-it,^)n+iO,0)) , (7) 

X 

where the summation over x corresponds to summing over all the spatial lattice sites, thereby 
projecting onto the momentum p = state. A 7r"*"7r"^ correlation function that projects onto 
the s-wave state in the continuum limit is 

a+.+ (p,t) = E Ee^'"^''''^^"(^'^)^"(^'y)^^(0,O)vr+(0,O)) , (8) 

|p|=p x,y 

where, in eqs. (jTj) and (jH)), vr+(t,x) = u(t, x)75(i(t, x) is an interpolating field for the tt"*". 

Steps were taken to optimize the overlap between the interpolating fields and the one- and 
two-pion hadronic states. First, the propagators calculated by LHPC, which all have sources 
centered about x = 0, were smeared [31, 40] in the neighborhood of x = to maximize 
overlap with the single-hadron states. Therefore, in eq. ((Tj) and eq. (jH)), 7r+(0,0) denotes 
an interpolating field constructed from smeared-source light-quark propagators. Second, 
we projected onto two-pion states that are perturbatively close to the energy eigenstates 
of interest. The two-pion states with zero total momentum that transform in the Ai rep- 
resentation of the cubic group are, in general, linear combinations of the non-interacting 
finite- volume eigenstates, 

|7r+7r+) = do |7r+(p = 0)7r+(p = 0)) + rfi E k+(p)^+(-p)) + - (9) 

where do and di are complex coefficients. In the absence of interactions, the ground state is 
given by the first term in eq. the first excited state by the second term, and so on. In 



^ Note that the source of each propagator was placed on the t = 10 time shce, which led to a maximum of 
At — 22 usable time slices. 
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the regime \a2\ <C I/, the effect of the interaction is small, as is clear from eq. and each 
of the terms in eq. (jH)) are approximate eigenstates. The momentum projection in eq. (jH)) 
makes it significantly easier to extract the energies of the interacting eigenstates. 

In order to determine the amplitudes and energies of the eigenstates that our correlation 
functions contain, we used a variety of fitting methods. These include standard covari- 
ant fitting, non-covariant fitting, the effective-mass method, singular-value decomposition, 
single-exponential fits, and multiple-exponential fits. Table IH] provides a summary of the 
results of our calculation. For each light quark mass we have chosen to extract the neces- 
sary quantities by fitting the correlation functions over the fitting interval given in Table ITTl 
Using a covariance-matrix, fitting procedure, along with Jackknifing over the lattice con- 
figurations, we obtain the central values and statistical errors (the first error quoted for each 
quantity) shown in Table |nj Varying the fitting interval over reasonable ranges provides an 
estimate of one of the systematic errors associated with the determination of these quanti- 
ties. The range of central values from a small number of different fitting intervals is used 
to determine the second error associated with m^r, /tt, t^-k/ f-K^ ^Eq, 02 , m^a2 and AEi, 
given in Table ITTl In the calculation of the remaining quantities in Table ITTl their associated 
systematic error was estimated by combining statistical and systematic errors of the input 
quantities in quadrature. For m^r, and |p| given in units of MeV, there is an additional 
systematic error arising from the scale setting, which we estimate from the uncertainty in 
the counterterm, I4, contributing to which is discussed subsequently. 



77 ° 
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FIG. 1: Log plots of the ratio of correlation functions, Gt^, defined in eq. (fT?!]) . The left panel 
shows the correlation function G7r7r(0, t) for the three sets of MILC configurations, each of which is 
dominated by the tttt ground-state. The right panel shows the logarithm of the absolute value of 
Gt^tj {2tt / L^t) for the heaviest quark-mass with and without subtraction of the (small) ground-state 
component. 



The pion mass was determined by fitting a single exponential to the pion correlation 
function, eq. (|7j). The point-smeared correlator (point-sink and smeared-source) was found 
to be statistically superior to the smeared-smeared correlator (where the sink is smeared with 
the same smearing function as the source) both for the one- and two-pion correlators. The 
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smeared-smeared correlator was used only in the determination of the pion decay constant, 
(see below). 

In the relatively large lattice volumes that we are using, the energy difference between 
the interacting and non-interacting two-pion states is a small fraction of the total energy, 
which is dominated by the mass of the pions. In order to extract this energy difference we 
followed Ref. and formed the ratio of correlation functions, G.„.„{j),t), where 

G^^M - %77# - E A e-^^" * , (10) 

'^■kV') „=o 

and the arrow becomes an equality in the limit of an infinite number of gauge configura- 
tions. In GTr^{p,t), some of the fluctuations that contribute to both the one- and two-pion 
correlation functions cancel, thereby improving the quality of the extraction of the energy 
difference beyond what we are able to achieve from an analysis of the individual correlation 
functions. 

In the p = case, where we project perturbatively close to the ground state, the correlator 
G'7r7r(0, t), shown in Fig.Hl can be fit by a single exponential beyond the first few time slices, 
and the ground-state energy difference, AEq, can be determined quite cleanly. For the 
first excited level, with p = 27r/L, the momentum projection in eq. (jH)) eliminates all but 
approximately 10% (in amplitude) of the ground state contribution, which contributes with 
opposite sign, as shown in Fig. ^ We analyzed this correlator in two ways. First, by 
fitting a single exponential to what remains after subtracting the ground-state contribution 
determined at a distant time slice. Second, by fitting two exponentials, either keeping the 
value of the ground state energy fixed (three parameter fit) or letting that vary also (four 
parameter fit). The difference between these procedures, as well as the difference coming 
from the somewhat arbitrary choice of fitting ranges, is significant only in the value of the 
two-pion energies and is incorporated in the estimated systematic error. The results of 
our analysis and the quantities relevant to the determination of the scattering length are 
summarized in Table ITTl 



IV. ANALYSIS 

A main focus of the analysis of lattice data is the extrapolation of the lattice values of var- 
ious physical and unphysical parameters that hadronic quantities computed on the lattice 
depend on, in order to make direct comparison with experimental hadronic quantities. Ex- 
trapolations in the light-quark masses, the finite lattice spacing and the lattice volume are 
currently required. For small enough quark masses and lattice spacings, and large enough 
volumes, one can rigorously perform such extrapolations using low-energy effective field the- 
ories. What defines small enough are the dimensionless quantities bA^ -C 1, m.„/A^ <^ 1 
and m^L 3> 1. The low-energy effective field theory that can be used to describe the vrvr 
scattering length is Y-P T supplemented to include the finite lattice spacing, based on the 
Symanzik action (41. 4^ . 



Our calculations have generated pions with masses that are less than ~ 500 MeV, which is 
likely an upper limit to the range of applicability of the chiral expansion. We are confident 
that the results of the calculations at the lightest two pion masses fall within the chiral 
regime and can be analyzed with X"PT, while the results of the largest mass point is at the 
edge of applicability of x-PT. 
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The effect of the finite lattice spacing is not universal and depends upon the particular 
discretization used. In our mixed-action calculation, both the sea quarks and valence quarks 
have been included with actions that have lattice spacing corrections beginning at (9(6^). 
While rrires 7^ implies 0{b) corrections, these effects are exponentially suppressed with 
the extent of the fifth dimension and are therefore safely ignored given our current level 
of precision. In addition, the sea-quark discretization effects are further suppressed by the 
strong coupling constant since the sea-quark action is perturbatively improved. A gener- 
alization of x-PT suitable to describe this combination of mixed-quark-action discretization 



has recently appeared in Ref. [43]. Our calculation has been performed at a single value 
of the lattice spacing and as such we are unable to systematically explore the impact of 
discretization errors on the determination of the scattering amplitude. In the analysis that 
follows, we have assumed that the lattice-spacing effects are much smaller than other sources 
of uncertainty in this calculation, and are neglected ^. 

While the power-law finite-volume dependence of the energy of the two-pion states is 
used to extract the scattering amplitude using Liischer's method, there are also exponentially 
suppressed, finite- volume contributions to the observables resulting from pions "going around 
the world". These corrections are proportional to e"'"'"^ and are further suppressed by a 
factor m^/A^, since they can arise only from pion loop diagrams. As each of the pion masses 
that we compute in this work satisfy rriT^L ^ 1, these exponentially-suppressed finite- volume 
contributions are neglected in the analysis. 



A. Scale Setting 

In order to determine quantities in physical units, the lattice scale must be traded for 
a physical scale. In many calculations the mass of the p-meson is used. While this seems 
fairly convenient ,the light-quark mass dependence of the p-meson is controversial (for recent 
work, see Ref. |^), and moreover, for light enough pions, the p-meson becomes unstable, 
and therefore determining its mass in the chiral regime requires a series of calculations of the 
energy levels of two-pion states in a finite- volume. A short-distance length scale is also com- 
monly used 0], the Sommer scale, tq. Such a short- distance scale depends logarithmically 
upon the light-quark masses, as stressed in Ref. 4^|, however lattice calculations indicate 



that the numerical size of the light-quark mass dependence is quite small ^49.,^], i.e. the co- 
efficients of unknown operators in the quarkonium chiral lagrangian are numerically small ^. 
In this work the pion decay constant, is used to set the physical length scale. This has 
the advantage that: (i) is a low-energy parameter with a well-known dependence on the 
quark masses (as well as the lattice volume and lattice spacing) determined by x-PT; (ii) 
the decay constant is easily computed on the lattice with domain-wall valence quarks. We 
calculate from pion correlation functions via |5ll | 



Asv_ 2 ^/2 ( rudy^f + rrires ) 
^55 ml 



r " •■OH — "aiii/ I ' ■ "res j / 1 1 \ 



^ In the / = TTTT channel, unlike the / = 2 tttt channel, diagrams involving the double pole in the ij' 
propagator, proportional to b^, lead to modifications to Liischer's formula in cq. 10} that are enhanced by 



powers of L |4J, |45 

^ The coeSicients are expected to scale as R^, where R is the characteristic size of the QQ system, from 
naive dimensional arguments. 
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where As-p is the amphtude of the one-pion correlation function resulting from a smeared- 
source and point-sink, and Ass is the amplitude which results from a smeared-source and 
smeared-sink. The values of m^,„ f and rrires used in this work can be found in Table HI 
The chiral expansion of /^r is j52] 



U = f 



1 + 



m 



I4 + 0{m^ 



(12) 



where / is the chiral-limit value of the decay constant, m is the pion mass at leading order 
in the chiral expansion; 7j = logA^/m^ with Aj an intrinsic scale that is not determined by 
chiral symmetry. In what follows we will denote the physical values of the various parameters 
with a p/iy-superscript. One can use eq. fll2|) to find 




fp 



phy 



1 + 



47r2 




(13) 



m 



phy'' 



phy 



if' + O 



ml 



as well as other variants of this formula that differ only by terms of higher order in the chiral 
expansion. By determining rriT^/ f^^ in a lattice calculation, and using the experimental values 
of rrtTT and /tt (and I4), the pion-decay constant in the lattice calculation at the value of the 
pion mass, fni^n/ fw), is determined at a given order (in this case next-to-leading order) 
in the chiral expansion. Higher-order corrections to this scale setting can be determined 
systematically in the chiral expansion. We use /^'^^ = 132 MeV and 771'^^^ = 138 MeV. The 

uncertainty in the scale-independent parameter ifl^ = 4.4 ±0.2 js^E^ introduces a source 
of systematic error of less than 2% in both the pion mass and decay constant. The lattice 
spacings obtained using this method are given in Table |n] and may be compared with the 
determination by MILC [i^, b = 0.1243 ± 0.0015 fm, using the Sommer scale. The values of 
the pion masses and decay constants in physical units computed in this paper are consistent 
with the same quantities computed on the same lattices by MILC in Ref. j49| . For instance, 
MILC finds = 299.99 ± 0.32 ± 3.4, 356.46 ± 0.27 ± 4.0, 494.27 ± 0.25 ± 5.6 on the 0.007, 
0.010 and 0.020 lattices, respectively, where the first error is statistical and the second arises 
from scale setting. One should keep in mind that (i) results differ at (9(6^) and 0{m^)] (ii) 
the tuning of the domain-wall fermion mass to reproduce the lightest staggered pion mass is 
done to order 1% accuracy; (iii) Our scale setting procedure includes electromagnetic effects 
via fP^y while MILC's does not. 

Rather than using the low-energy constant if' as an input to the scale setting procedure, 
one can perform a two-parameter fit of the lattice data in Table HTl to the form 





iphy 

''4. fit 



(14) 



where fj'^' is the pion decay constant in lattice units. This procedure gives bfn 

-phy 



0.1274 ± 0.0007 ± 0.0003 fm and / 



4 fit 



4.412 ± 0.077 ± 0.068 where the first error is 
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statistical and the second is an estimate of tlie systematic error. Tlius tliis scale set- 
ting procedure is remarkably robust and consistent. One may wonder about the rel- 
evance of the chiral logarithm. Repeating the fitting procedure with /J'"' {mT^/ f^^) = 



'Jfit 



jphy 



■fit 



yields b 



'fit 



0.1330 ± 0.0001 ± 



0.0001 fm and Cfit = 1.407±0.010±0.009, which are not consistent with MILC scale setting 



or the experimental value of l^^^ , respectively. It would appear that the chiral logarithm is 
resolved by our data at this order in the chiral expansion. 



B. The Scattering Length 

With small quark masses and momenta, tttt scattering can be reliably computed in X"PT. 
The leading-order result (equivalent to current algebra) was computed in Ref . [s^l , and the 
one-loop TTTT amplitude was computed in Ref. While this amplitude is now known at 
the two-loop level jH, , given our current lattice data, we choose to analyze our lattice 
results at one-loop level. The one-loop expression for the J = 2 vrvr scattering length is 



2 

mi 



2 



(15) 



where Itt-k{ij) is a linear combination of scale-dependent low-energy constants that appear 
in the 0{p^) chiral lagrangian [H^ (see Appendix We define I-,,-,, = /,r7r(/^ = 47r/^), 
and therefore we can simply use the ratio m^/ f^, computed on the lattice to determine the 
scattering length using eq. (fT3j) . The difference between using the lattice and a fixed 
in the argument of the logarithm modifies the scattering length only at higher orders in the 
chiral expansion. 

The lowest-lying energy eigenvalues in the lattice volume, shown in Table |nj allow us to 
determine the / = 2 tttt scattering lengths at the different light-quark masses via eq. (jH). Our 
results for the scattering lengths, and other parameters are presented in the summary table, 
Table ITTl The location of the first excited state in the lattice volume allows, in general, for a 
determination of the phase-shift at non-zero values of the pion momentum via eq. (^. For 
the lattice parameters in these calculations we were able to extract the 1 = 2 tctc phase-shift 
at one (large) momentum at the largest quark mass, which is shown in Table |n] For the two 
lighter quark masses, the first excited state is very near the four-pion inelastic threshold, 
and a simple extraction of the tttt phase-shift is not possible. 

The results of our calculation of the product 7717^02 are shown as a function of m^/ f^^ in 
Fig. Q. In addition, we have shown the lowest pion mass datum from the dynamical calcu- 
lations of the CP-PACS collaboration 2^1 ^. The uncertainty in the CP-PACS measurement 



is significantly smaller than that of our calculation and the agreement is very encouraging. 
In order to extrapolate mT,a2 to the physical value of 171^^/ f.^-, we performed a weighted fit of 
eq. (fT3j) to the three data points in Table ITTl and extracted a value of the counterterm Z^^r. 
As both quantities, 7717^02 and m^//^, are dimensionless there is no systematic uncertainty 
arising from the scale setting (I4). We determined that /tttt = 3.3 ± 0.6 ± 0.3, where the first 



^ We have shown the CP-PACS data point at the hghtest pion mass, and at the smallest lattice spacing, 
P = 2.10, and have not attempted to extrapolate their result to the continuum. This lattice spacing is 
comparable to the one used in this work. 
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TABLE II: The summary table. The central value for each quantity is determined by fitting the 
appropriate correlation function over the indicated "Fit Range" . The first uncertainty is statistical. 
The second is an estimate of the systematic error in the fitting process (including varying the fitting 
range). For quantities with units of MeV, the third uncertainty is theoretical, and is due to the 
uncertainty in the low-energy constant, I4 , originating from scale-setting. The chi-square test of 
fit refers to the extraction of AEq. 



Quantity 


mi = 0.007 


mi = 0.010 


mi = 0.020 


Fit Range 


6-13 


5-15 


7-15 


(1-u.) 
(l.u.) 


n 1 onn -i- n nnoi -i- n nno 
u.iyuu It u.uuzi zt u.uuz 

0.0937 ± 0.0012 ± 0.001 

2.030 ± 0.040 ± 0.03 


n 99/1 *^ -1- n nni n -i- n nnn^ 

U.ZZ^to It U.UUiU It u.uuuo 

0.0959 ± 0.0007 ± 0.0002 
2.338 ± 0.022 ± 0.005 


n Qi Qi -1- n nni 9 -i- n nni 7 

U.OiOi It U.UUiZ It U.Uui ( 

0.1021 ±0.0007 ±0.0012 
3.065 ± 0.024 ± 0.030 


X /d.o.f. 


0.19 


0.84 


1.03 


AEo (l.u.) 
a2 (l.u.) 

mT,a2 
C 


0.0109 ± 0.0013 ± 0.0003 

-1.12 ±0.12 ±0.02 
-0.212 ±0.024 ±0.004 
0.29 ±0.16 ±0.05 


0.0080 ± 0.0005 ± 0.0003 

-0.99 ±0.06 ±0.04 
-0.222 ±0.014 ±0.009 
0.021 ±0.067 ±0.041 


0.0073 ± 0.0007 ± 0.0004 
-1.22 ±0.09 ±0.07 
-0.38 ±0.03 ±0.02 
0.017 ±0.082 ±0.057 


U (MeV) 
m„ (MeV) 
b (10-2fm) 


144.7 ± 0.5 ± 0.4 ± 1.0 
293.7 ±5.9 ±4.4 ±2.0 
12.78 ±0.19 ±0.18 ±0.09 


148.8 ±0.3 ±0.1 ±1.5 

347.9 ±3.3 ±0.8 ±3.4 
12.72 ±0.11 ±0.09 ±0.13 


158.0 ±0.3 ±0.4 ±2.8 
484.4 ±3.9 ±4.9 ±8.5 
12.75 ±0.13 ±0.21 ±0.23 


A^i (l.u.) 

IpI (MeV) 
5{p) (degrees) 


0.482 ± 0.070 ± 0.049 
607 ± 44 ± 42 ± 15 


0.390 ± 0.030 ± 0.035 
551 ±29 ±31 ± 12 


0.308 ± 0.009 ± 0.005 

544 ± 14 ± 12 ± 10 
-43 ± 10 ± 5 



error is statistical and the second is an estimate of the systematic error (see Fig. Q). This 
fit of /tttt allows, through eq. (|T5|l . a prediction of the scattering length at the physical value 
of the light-quark masses, which we find to be 

m^as = -0.0426 ±0.0006 ±0.0003 ±0.0018 . (16) 

The last uncertainty, ±0.0018, is the largest and is an estimate of the systematic error 
resulting from truncation of the chiral expansion of the scattering length ®. The two-loop 
expression for the scattering length jssl . Isq is given by 

.2 



ml 



1 ± 



3mi 



167r2/2 



± 



647r4/.^ 



TTl 

log^ ± i^M 



, (17) 



where l^} and l^^J are linear combinations of undetermined constants that appear in the 
C(p^) and 0{p^) chiral lagrangians 5^, 5^ (see Appendix IXJ . It is not possible to provide 



As a consistency check, one may fit the leading-order x-PT result multiplied by an arbitrary coefficient: 

2 

m-7ra2 = —A -^ji- This yields A = 1.058 ± 0.059 ± 0.032 where the first error is statistical and the second 
is an estimate of the systematic error. The scattering length at the physical value of the light-quark 
masses is, in this case, m^oa = -0.0462 ± 0.0026 ± 0.0014 ± 0.0018. 
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m I f 

FIG. 2: The results of this lattice QCD calculation of m-,^a2 as a function of rriTr/ fn (ovals) with 
statistical (dark bars) and systematic (light bars) uncertainties. Also shown are the experimental 
value from Ref. [s^l (diamond) and the lowest quark mass result of the dynamical calculation of 
CP-PACS (square). The gray band corresponds to a weighted fit to our three data points 
using the one-loop x-PT formula in eq. (fT3|) which gives /tttt = 3.3 ± 0.6 ± 0.3 (the shaded region 
corresponds only to the statistical error). The dashed line is the tree- level X"PT result. 



a meaningful fit of these three undetermined constants from the few data points we have 
calculated. While there are estimates of these low-energy constants from a variety of sources, 
using these estimates in our extrapolation to the physical point would amount to trading 
one unknown systematic error for another. In order to estimate the systematic error due 
to extrapolation, we first set /^^^ = ifl = and refit /tttt keeping only the double-log piece 
at two-loop order. We then set l^^} = and do a two-parameter fit to 1^,^^ and l^J. The 
difference between the extrapolation of the one-loop expression and the extrapolation of 
the two-loop expression with these simplifications gives the estimate of the systematic error 
resulting from truncating the chiral expansion. 

In Fig. (j2I) we show the unique prediction of leading order X"PT for the scattering length 
(the dashed line), which is seen to agree remarkably well with both the lattice data and the 
physical value. The gray band in Fig. Q shows the statistical la region. 

In order to better isolate the contributions from higher orders in X"PT it is convenient 
to define a scale-independent "curvature" function C: 

Mfia, , _ 3m; ( ^ ^\ _ (18, 
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where, by construction, C = at tree level. The values for C(m7r//7r, ^tttt) that we have 
calculated are listed in Table ITTl and are shown in Fig. Q as a function of ra-^j f.^ together 
with the experimental point and the CP-PACS result. A weighted fit of eq. ()18p to the 



0.5 



0.4 



0.3 



0.2 



0.1 



-0.1 



-0.2 



^ Experiment 

■ CP-PACS 
• This Work 



m If 

FIG. 3: C{mT^/ f-,^^ /tttt) as a function of m-,^/ for our three lattice data points (ovals) with statistical 
(dark bars) and systematic (light bars) uncertainties, plotted with the experimental value from 
Ref. ^3] (diamond) and the lowest-mass dynamical CP-PACS [i^l point (square). The gray band 
corresponds to a fit to our three data points (weighted by statistical errors) using the one-loop 
X-PT formula in eq. (|18j) . 



results of the lattice calculation gives: /^,r = 3.3 ± 0.6 ± 0.2, consistent with the fit to m^a2. 
The gray band in Fig. (jS)) shows the la region for C. The extrapolated value of the scattering 
length is obviously the same as in the direct fit to 771,^(32 • 

The value of /,r7r favored by the fits are such that there is an almost perfect cancellation 
between the counterterm and the logarithm in eq. (|15|) in the range of m^r considered. This 
cancellation may be an unfortunate coincidence in this channel. A more refined lattice QCD 
calculation is required in order to detect the predicted chiral curvature. 

The best experimental determination of the 1 = 2 scattering length is obtained through 
an analysis of K{e4) decays HJ, which gives m^az = -0.0454 ± 0.0031 ± 0.0010 ± 0.0008 
where the first error is statistical, the second is systematic and the third is theoretical. This 
data point is plotted in Fig. Q and in Fig. Q. The two- loop x-PT "prediction" [s^l is 
= —0.0444 ± 0.0010. Our result is in very good agreement with these determinations. 
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C. The Phase Shift 



We are only able to extract the phase-shift at one non-zero value of the pion momentum, 
and only at the heaviest pion mass, despite having a relatively clean signal for the first 
excited state in the lattice volume for all three sets of MILC configurations. The reason for 
this is that the pion masses are sufficiently light that the first excited state is at an energy 
very near the four-pion inelastic threshold on the two sets of configurations with the lightest 
pion masses. At m^r ~ 484 MeV, and for pion momentum |p| ~ 544 MeV, the phase-shift is 
found to be 5 = —43 ± 10 ± 5 degrees. 

V. CONCLUSIONS 

In this paper we have presented the results of a lattice QCD calculation of the / = 2 tttt 
scattering length performed with domain-wall valence quarks on asqtad-improved MILC 
configurations with 2+1 dynamical staggered quarks. The calculations were performed at 
a single lattice spacing of 6 ~ 0.127 fm and at a single lattice size of L ~ 2.5 fm with 
three values of the light quark masses, corresponding to pion masses of m^r ~ 294, 348, and 
484 MeV. We have also presented the phase-shift at the heavier quark mass. 

We have used one-loop X"PT to fit the combination of counterterms contributing to vrvr 
scattering at next-to-leading order to the lattice data and extrapolated in the light-quark 
masses down to the physical point. At the one- loop level we are able to make a prediction 
for the value of the scattering length, m^a2 = -0.0426 ± 0.0006 ± 0.0003 ± 0.0018, which 
agrees within errors with the experimental value. 
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APPENDIX A: THE 1 = 2 SCATTERING LENGTH AT TWO LOOPS 



To accuracy 0{m^) in X"PT (two loops), the scattering length is given by j& 



where 
B 



2 + 



mz 



262 227r2 ^ „ 

+ 647r2 i3 

9 9 



(Al) 



3 m 3 

T log ^ + 8(li+ h) + 2(l3- h) ^ 



+2 



/2 Il67r2 



47 1 1861 

24 /i + 16 /2 + 15 /3 - 6 /4 + ■ TTT^ log ^ + 



31 



6 \ 167r2 



1 TTi^ \ 1 TTl 

log^ - — [/4 + 2(2/i + Z3 + 6/2)]l0g^ 



12 16772 ^ /i2 23047r2 

2 
7r 
2" 



-5/^ - 8ll + 4/4 [8/1 + 3/3 + 8/2] + ri + 16r4 | 



(A2) 



The h and are scale -dep endent low-energy constants that appear in the 0{p'^) and C(p^) 
chiral lagrangians 5^, IJ, respectively. 
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